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ABSTRACT 

We report on spectral and timing observations of the nearest millisecond pulsar 10437^715 with the Chandra 
X-ray Observatory. The pulsar spectrum, detected up to 7 keV, cannot be described by a simple one-component 
model. We suggest that it consists of two components, a nonthermal power-law spectrum generated in the pulsar 
magnetosphere, with a photon index 7 « 2 , and a thermal spectrum emitted by heated polar caps, with a tempera¬ 
ture decreasing outwards from 2 MK to 0.5 MK. The lack of spectral features in the thermal component suggests 
that the neutron star surface is covered by a hydrogen (or helium) atmosphere. The timing analysis shows one 
X-ray pulse per period, with a pulsed fraction of about 40% and the peak at the same pulse phase as the radio 
peak. No synchrotron pulsar-wind nebula is seen in X-rays. 

Subject headings: stars; neutron — pulsars: individual (PSR 10437^715) — X-rays: stars 


1. X-RAYS FROM MILLISECOND PULSARS: THERMAL VERSUS 
NONTHERMAL 

Millisecond (recycled) radio pulsars are distinguished from 
ordinary pulsars by their very short and stable periods, P ^ 

10 ms, P ~ 10“^' - 10“*® s s“*. It is generally accepted that 
they are very old objects, with spin-down ages r = P/2P 
10®- 10*** yr and low surface magnetic fields B oc (PP)^!'^ 
10^-10*® G (e.g., Taylor, Manchester, & Lyne 1993). Similar 
to ordinary pulsars, a millisecond pulsar can emit nonthermal 
X-rays from its magnetosphere, with a hard power-law spec¬ 
trum and sharp pulsations. In addition to this nonthermal ra¬ 
diation, thermal X-rays can be emitted from the neutron star 
(NS) surface, provided the surface is hot enough. According to 
the models of NS thermal evolution (albeit rather uncertain at 
these old ages), recycled pulsars are too cold (surface tempera¬ 
ture r ^ 0.1 MK — see, e. g., Tsuruta 1998) to be detectable 
in X-rays. However, their polar caps can be heated up to X-ray 
temperatures by relativistic particles impinging onto the mag¬ 
netic poles from the acceleration zones in the magnetosphere. 
The radio pulsar models (e.g., Cheng & Ruderman 1980; Arons 
1981; Michel 1991; Beskin, Gurevich, & Istomin 1993) predict 
polar cap radii Rpc ~ /Pc)^/^ (where R « 10 km is the NS 

radius), i.e., Rpc ~ 1-5 km for millisecond pulsars, although dif¬ 
ferent models predict quite different polar cap temperatures, in 
the range of 1-10 MK. Detection of the polar cap thermal radi¬ 
ation would allow one to discriminate between various models 
of radio pulsars, study the properties of NS surface layers, and 
constrain the NS mass-to-radius ratio (Pavlov & Zavlin 1997; 
Zavlin & Pavlov 1998 [ZP98]). However, just as in the case 
of ordinary pulsars, this radiation is detectable only if it is not 
buried under stronger nonthermal radiation. The current the¬ 
oretical models are not elaborate enough to predict in which 
(if any) of millisecond pulsars the thermal component can be 
brighter than the nonthermal one (in particular, both the ther¬ 
mal and nonthermal luminosities are expected to increase with 
spin-down energy loss E, perhaps with different rates). There¬ 
fore, we have to rely upon the analysis of X-ray observations to 
distinguish the thermal and nonthermal components. 


The X-ray observatories ROSAT, ASCA, and BeppoSAX 
have detected 11 millisecond pulsars (nearly 1/3 of all X-ray- 
detected rotation-powered pulsars — see Becker & Pavlov 2001 
for a recent review). Live of these pulsars are identified in 
X-rays only by positional coincidence with the radio pulsars 
and, due to the low number of recorded counts, provide only 
crude flux estimates. The radiation from 3 pulsars — B1821- 
24 (Saito et al. 1997), B1937H-21 (Takahashi et al. 2001), and 
J0218H-4232 (Mineo et al. 2000) — is clearly nonthermal: their 
power-law spectra, detected with ASCA and BeppoSAX up 
to energies of 5-10 keV, are very hard, with photon indices 
7 ~ 1, and their pulse profiles show sharp peaks. Interestingly, 
these 3 pulsars are characterized by particularly large E values, 
E = (2-20) X 10^^ erg s“*, and their magnetic fields at the light 
cylinder, Bic = B(R/Ric)^ ~ 10® G, are close to that of the Crab 
pulsar. The case for the other 3 pulsars —10437^715 (Becker 
& Trumper 1993, 1999 [BT93, BT99]; ZP98), J2124-3358 
(BT99), and J0030H-0451 (Becker et al. 2000) — is less cer¬ 
tain. These pulsars show broad peaks of X-ray pulsations, but 
it does not necessarily mean that their radiation is thermal be¬ 
cause broad peaks can be produced by nonthermal emission at 
some viewing angles. High-quality spectra have been recorded 
for the brightest of these pulsars, 10437^715, but their inter¬ 
pretation has been controversial — e.g., ZP98 suggest that the 
radiation detected with ROSAT and EUVE can be interpreted 
as thermal radiation from hot polar caps, whereas BT99 argue 
that the radiation is nonthermal (see §2). To resolve this con¬ 
troversy, the pulsar needed to be observed at energies above the 
soft ROSAT and EUVE bands {E'^2 keV), and with high spa¬ 
tial resolution to avoid contamination of the pulsar emission by 
a nearby AGN which compromised the ASCA and BeppoSAX 
data. 

The Chandra X-ray Observatory provides both the superb 
spatial resolution and high throughput at higher energies, to¬ 
gether with timing capability. In this paper we present the re¬ 
sults of our observations of PSR 10437^715 with Chandra. We 
start from a summary of the previous results on PSR J0437- 
4715 in §2. The spectral and timing analyses of the Chandra 
data are presented in §3 and §4. Implication of the results are 
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discussed in §5. 

2. PREVIOUS X-RAY OBSERVATIONS OF PSR 10437^715 

At a distance d = 139 ±3 pc (van Straten et al. 2001), 
PSR 10437^715 is the nearest and brightest X-ray millisec¬ 
ond pulsar known at both radio and X-ray wavelengths. It was 
discovered by Johnston et al. (1993) during the Parkes South¬ 
ern radio pulsar survey. It has a spin period P = 5.76 ms, and, 
after correcting the observed period derivative for the kinematic 
term due to its large proper motion, a characteristic age r ~ 4.9 
Gyr, magnetic field B ~ 3 x 10^ G, and rotation energy loss rate 
E ~ 3.8 X 10^^ erg s“'. It is in a 5.74 d binary orbit with a white 
dwarf companion of a low mass M ~ Q.2Mq. Optical obser¬ 
vations in Hq, have revealed a bow-shock pulsar-wind nebula 
(PWN), caused by the supersonic motion of the pulsar through 
the interstellar medium, with the bow-shock apex at about 10 " 
south-east of the pulsar, in the direction of the pulsar’s proper 
motion (Bell, Bailes, & Bessell 1993; Bell et al. 1995). 

PSR 10437^715 was the first millisecond pulsar detected in 
X-rays — BT93 observed it with the ROSAT Position Sensi¬ 
tive Proportional Counter (PSPC) in 1992 and discovered X- 
ray pulsations with a single broad pulse and a pulsed fraction 
/p ~ 30% in the 0.1-2.4 keV energy range. They also found, 
and Halpern, Martin, & Marshall (1996; HMM96 hereafter) 
later confirmed, that /p apparently increases with photon en¬ 
ergy E and reaches about 50% at E = 0.6- 1.1 keV. However, 
BT99 found no energy dependence of /p from the PSPC obser¬ 
vation of 1994. On the other hand, ZP98 analyzed the same data 
and claimed that the energy dependence of the pulse profile is 
statistically significant. 

The pulsar spectrum is incompatible with a simple blackbody 
(BB) model even in the narrow PSPC band (BT93), but it can 
be fitted with a power-law (PL) model, with a photon index 
7 = 2.2-2.5 (BT93; HMM96; ZP98; BT99), indicating that the 
soft X-ray radiation might be nonthermal. However, the proba¬ 
ble energy dependence of the pulsed fraction makes this inter¬ 
pretation questionable. In addition, the combined analysis of 
the PSPC and EUVE Deep Survey Instrument (DSI) data show 
that the PL fit is only marginally acceptable (HMM96; ZP98). 

BT93 suggested that the pulsar spectrum consists of two 
components — in addition to a PL (7 Ri 2.9), representing 
magnetospheric emission, there is a BB component with Tbb ^ 
1.7 MK emitted from a small area of 0.05 km^, implying the 
existence of polar cap(s) on the NS surface. An alternative two- 
component model was explored by ZP98. They suggest that 
the temperature distribution in a polar cap should be nonuni¬ 
form because the heat, released by decelerating magnetospheric 
particles in subphotospheric layers, propagates along the sur¬ 
face out from the small polar area where the energy is de¬ 
posited, which results in a larger heated area with the tem¬ 
perature decreasing outwards. Applying the hydrogen NS at¬ 
mosphere models of Zavlin, Pavlov, & Shibanov (1996), ZP98 
show that the entire soft X-ray spectrum and the pulse profiles 
detected with ROSAT and EUVE can be interpreted as thermal 
emission from two opposite polar caps where the temperature 
drops from Tore =1-2 MK in a central cap (“core”) of radius 
Rcore = 0.2-0.5 km to Trini = 0.3-0.5 MK in the annulus (“rim”) 
with an outer radius Brin, = 2-5 km. The bolometric luminosity 
of the polar caps^, Lboi = (2- 3) x 10^° erg s“', is consistent with 


the prediction of the slot-gap model of radio pulsars developed 
by Arons (1981). The hydrogen column density towards the 
pulsar, hh.io = «h/( 10^° cm“^) = 0.1 -0.3, obtained by ZP98, 
is consistent with the properties of the interstellar medium in¬ 
ferred from observations of other objects in the vicinity of the 
pulsar. 

The ROSAT energy band, 0.1-2.4 keV, is too soft to discrim¬ 
inate between different two-component models — e.g., BT99 
have shown that equally good fits of the ROSAT PSPC spec¬ 
trum can be obtained with a two-component BBh-BB model 
or a broken PL. The pulsar has been observed at higher ener¬ 
gies with ASCA. Kawai, Tamura, & Saito (1998) report that 
the ASCA spectrum is more consistent with a BB model with 
7bb ~ 3 MK than with a PL model. However, because of poor 
spatial resolution of ASCA, the pulsar spectrum was heav¬ 
ily contaminated by the radiation from a neighboring Seyfert 
galaxy, that hampered the spectral analysis. Therefore, to in¬ 
vestigate the pulsar spectrum at higher energies and measure 
its pulse profile with high accuracy, we proposed to observe 
PSR 10437^715 with Chandra. The results of the spectral, 
timing, and spatial analyses are presented below. 


3. ACIS OBSERVATION AND SPECTRAL ANALYSIS 

PSR 10437^715 was observed with the Advanced CCD 
Imaging Spectrometer (ACIS; Garmire et al. 1992; Bautz et 
al. 1998) on 2000 May 29-30, with a 25.7 ks effective exposure. 
It was imaged on the back-illuminated chip S3 of the ACIS 
spectroscopic array. To mitigate the pile-up effect®, the pul¬ 
sar was offset by 3(9 from the optical axis and observed using a 
1/8 subarray, with a 0.4 s frame time. We extracted source-plus- 
background counts from a 5"-radius circle centered on the pul¬ 
sar image, which includes about 98% of pulsar counts. To mea¬ 
sure the background, we extracted counts from several regions 
at distances of at least 10" from the pulsar, where the pulsar 
contribution is negligible. The background, 0.020 counts ks“* 
arcsec”^, is rather uniform throughout the observed 8' x 1' sky 
area — in particular, we see no structures which could be in¬ 
terpreted as a PWN associated with the bow shock observed in 
Hq. The source count rate is 306 ±4 counts ks“' (0.12 counts 
per frame), for the standard grades 02346, in the 0.2-7.0 keV 
range (where the source prevails over the background). 

To examine the effect of pile-up on the pulsar spectrum, we 
extracted the count-rate spectra from a 1 "-radius circle centered 
at the pulsar (which contains 56% of source counts) and from 
an annulus with 1" and 5" radii. We compared the spectral 
shapes using a test and found them to be the same with a 
probability of 92%-95% (depending on binning chosen). Since 
the annulus spectrum should not be affected by pile-up, we con¬ 
clude that the effect of pile-up on the whole pulsar spectrum is 
negligible. 

To search for possible spectral features in the ACIS spec¬ 
trum, we used a method applied by Pavlov et al. (2001) to the 
dispersed X-ray spectrum of the Vela pulsar, obtained with the 
Chandra Low Energy Transmission Grating. We first binned 
the original detector channels of a 5 eV width in bins of an 80 
eV width (comparable with the intrinsic detector resolution of 
^100 eV). Then, we combined 4-5 sequential binned channels 
in groups of 320^00 eV widths to estimate the deviation in 


^The values of polar cap radii and luminosities given in ZP98 should be multiplied by and tt, respectively, because they used an erroneous normalization 
factor. This error does not affect their qualitative conclusions. 

^See Chandra Proposers’ Observatory Guide (POG); Tttp: / / asc. harvard. edu/udocs/docs/docs . html 
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the number of source counts in each binned channel of a given 
group from the mean value in the group. The maximum devia¬ 
tion does not exceed a 1.8 ct level, which allows us to conclude 
that the ACIS spectrum of PSR J0437^715 reveals no signifi¬ 
cant spectral features. For further analysis, we binned the pulsar 
spectrum in 75 energy bins in the 0.2-7.0 keV range. For fit¬ 
ting the ACIS data, we used CCD responses generated with the 
CIAO 2.1 software and CALDB 2.7 calibration files. 

To investigate as broad energy range as possible, we analyzed 
the ACIS spectrum together with the ROSAT PSPC spectra ob¬ 
served on 1992 September 20-21 (5.9 ks exposure) and 1994 
July 2-6 (9.9 ks exposure). (A detailed description of the PSPC 
data and their analysis can be found in ZP98.) Our attempts 
to simultaneously fit the ACIS and PSPC spectra immediately 
showed that they disagree with each other at energies below 
0.6 keV — models which fit the ACIS and PSPC spectra at 
E > 0.6 keV are well above the observed ACIS data points at 
lower energies (particularly in a 0.3-0.5 keV range). We have 
encountered the same problem with the Vela pulsar (Pavlov et 
al. 2001) and PSR 1055-52 (Sanwal et al. 2000) — the ACIS 
and PSPC spectra of these pulsars are in good agreement only 
if we discard ACIS events below ~ 0.6 keV. This inconsistency 
is apparently associated with errors in the ACIS response at 
low energies, so we have to use the ACIS data only in the 0.6- 
7.0 keV range (63% of the source counts in 50 energy bins). 

3.1. Single power law fit 

Fitting the separate ACIS spectrum in the 0.6-7 keV range 
with a PL model yields a best-fit photon index 7 = 4.1 (min 
Xy = 61.17/47 = 1.30), drastically different from that obtained 
for the PSPC spectrum (7 = 2.4). Not surprisingly, a single PL 
does not fit these spectra together (min = 366.8/96 = 3.82), 
clearly indicating that more components are needed to get a sat¬ 
isfactory combined fit. 


variations, apparently including polarization mode transitions 
(Navarro et al. 1997). However, it is reasonable to assume that 
the rotation axis of the pulsar is close the orbital momentum 
axis whose inclination angle, « 43°, was determined by van 
Straten et al. (2001). On the other hand, the radio pulse mor¬ 
phology and polarization strongly suggest that the magnetic in¬ 
clination a is very close to (/, probably the difference does not 
exceed 2°. We choose C = ct = 45° for our models of polar cap 
radiation. We have checked that variation of the angles by ±10° 
around these values makes no appreciable effect on the spectral 
parameters, and changes of the inferred radii of the polar caps 
do not exceed 15%. 

This model provides an acceptable fit to the PSPC and ACIS 
data (see Fig. 1, left panels), with a minimum x^ = 1-25 (94 
d.o.f) and best-fit parameters hh.zo = 0.6 ± 0 . 2 , T/c = 1.80 ± 
0.05 MK, Rpc = 0.15 ±0.02km (at c/ = 140 pc), and 7 = 2.8 ±0.1 
(the uncertainties here and below are given at a 68 % level; 
the polar cap parameters are given as measured at the NS sur¬ 
face; the gravitational effects are taken into account to calculate 
the model spectra as seen by a distant observer — see Zavlin, 
Shibanov, & Pavlov 1995). The bolometric luminosity of two 
polar caps is Lboi = (0.9 ± 0.2) x 10^*^ erg s“' ^ 0.2 x 10“^£, 
and the luminosity of the nonthermal component in the 0 . 1 - 
10 keV range is Lx = (3.1 ±0.4) x 10^® erg s“' 0.8 x 10“^£. 

The thermal component prevails in the 0.6-2.7 keV band; the 
PL component, dominating at higher and lower energies, de¬ 
termines the estimated value of uh and gives the main contri¬ 
bution to the EUVE DSI flux. Comparison with the EUVE 
DSI data shows that this model is only marginally acceptable, 
similar to the case of the single PL fit to the PSPC spectrum, 
described by HMM96 and ZP98. Also, the nu values allowed 
by the fit to the DSI, PSPC and ACIS data still exceed the val¬ 
ues nH ,20 = 0.1 - 0.3 derived from independent measurements of 
the hydrogen column density towards objects in the vicinity of 
PSR 10437^715 (see ZP98). 


3.2. Power law plus single-temperature polar caps 

It has been shown by ZP98 that even the narrow-band PSPC 
spectrum cannot be interpreted in terms of a single-temperature 
polar cap radiation — there is always an excess of observed 
counts over the model spectrum at higher energies, even for 
the hydrogen-atmosphere model that provides the hardest spec¬ 
trum. 


3.2.1. Hydrogen polar caps 

To explore the possibility that the observed X-ray emis¬ 
sion consists of thermal (polar cap) and nonthermal (magne- 
tospheric) components, we start with a two-component model; 
PL plus two polar caps covered with hydrogen or helium^. 
We assume the gravitational redshift parameter gy = (1 - 
2GM iRc^yl'^ = 0.769, which corresponds to the NS radius 
R= 10 km at M = 1 AMq. To model the polar cap component, 
we should specify some values for (/, the angle between the line 
of sight and rotation axis, and a, the angle between the mag¬ 
netic and rotation axes. Measuring these angles from the swing 
of the radio polarization position angle at the pulse center is 
difficult for this pulsar because of the complex position angle 


3.2.2. Blackbody polar caps 

If we assume a standard blackbody model for the polar cap 
radiation, we obtain qualitatively similar results, with a higher 
temperature of the thermal component, Lbb — 2.9 MK, and a 
much smaller radius of emitting area, Rbb — 0.04 km. Spec¬ 
tral fits with the simplified blackbody model consistently pro¬ 
duce higher temperatures and smaller sizes because X-ray spec¬ 
tra emerging from light-element atmospheres are harder than 
blackbody spectra (Zavlin et al. 1996). 

3.2.3. Iron polar caps 

Although the pulsar spectral fits are very comfortable with 
featureless models, it is worthwhile to check a hypothesis that 
the surface layers of polar caps are comprised of heavier chem¬ 
ical elements, e. g., iron. Using NS iron atmosphere models 
(Rajagopal & Romani 1996; Zavlin et al. 1996*) for the thermal 
(polar cap) component, we found that the single PSPC spec¬ 
trum fits well with a power law plus iron atmosphere model, 
whereas the ACIS data of much better spectral resolution are in¬ 
consistent with this model, particularly around the (redshifted) 
L and K absorption features in the iron atmosphere spectrum 


'The difference between the hydrogen and helium spectral models is very small in the X-ray range, at the temperatures of interest. 

^The iron atmosphere spectra presented by Zavlin et al. (1996) were computed with an error in interpolating the opacity tables. We are grateful to R. Romani who 
brought our attention to this issue. The error has been corrected in the models used in this paper. An example of a corrected model spectrum has been presented by 
Pavlov & Zavlin (2000). 
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(at energies of 0.85 keV and 4.9 keV for gr = 0.769 — see right 
panels in Fig. 1). Varying the gravitational redshift shifts the 
feature energies in the model spectrum, but the features remain 
too strong to be consistent with the smooth observed spectrum 
at any allowed values of gr. 

3.3. Power law plus polar caps with nonuniform temperature 

ZP98 showed that a polar cap model with a nonuniform tem¬ 
perature distribution along the cap’s surface fits the PSPC spec¬ 
trum without invoking a nonthermal component. They sug¬ 
gested a two-step approximation for the polar cap temperature 
distribution, a hot “core” plus a colder “rim” (see §2 and left 
panels of Fig. 2). Applying this model to the ACIS data shows 
that the observed spectrum significantly exceeds the model at 
E >2 keV (the upper limit of the PSPC energy range) and ob¬ 
viously requires one more component. A natural choice for the 
additional component is a PL. Adding the nonthermal compo¬ 
nent greatly improves the fit (right panels in Fig. 2) reducing 
xj, from 2.39 to 1.13. This model has seven fitting parame¬ 
ters, some of them are strongly correlated (ZP98). The cor¬ 
relation is particularly strong between hh and Rnm- Since the 
best-fit values of the hydrogen column density are well consis¬ 
tent with the values nn.io = 0.1-0.3 obtained from the indepen¬ 
dent estimates, we chose to fix the column density to weaken 
the correlations between the other six fitting parameters. For 
nH ,20 = 0.2, we obtain Tiore = 2 . 1 !o ;3 ^rim = 0.54+g;^^ MK, 
Rcoie = 0.km, Rrim = 2 . 0^0 2 km, and 7 = 2.2+°;^. Varying 
the hydrogen column density within the nH ,20 = 0.1 -0.3 range 
results in a minor change (< 7%) in the values of Rnm and Rcore 
(at the other parameters fixed). In this model the bolometric lu¬ 
minosity of two polar caps is Lboi = (2.3 ± 0.6) x 10^° erg s“' 
0.6 X lO"^^, and the PL luminosity in the 0.1-10 keV range is 
Lx = (0.7 ± 0.3) X 10^° erg s“' ~ 0.2 x 10“^£. Note that this 
model is in good agreement with the EUVE DSl data. 

One may assume that the additional component to fit the hard 
tail of the ACIS spectrum is also a thermal one, instead of the 
PL. However, this thermal component has to be of a very high 
temperature, T ~ 12- 15 MK, emitted from a very small area 
of 1-2 m radius. Although we are not aware of reliable cal¬ 
culations of temperature distributions around pulsar’s magnetic 
poles, we believe that such a steep gradient in the temperature 
distribution is unlikely. 

3.4. Broken power law 

Although the power-law plus polar-cap models fit the 
PSPCh-ACIS spectra satisfactorily, one cannot rule out a pri¬ 
ori that the (purely nonthermal) pulsar spectrum bends down 
somewhere in the X-ray band, which might be caused, e.g., 
by a deficit of high-energy electrons in the pulsar’s magne¬ 
tosphere. To test this hypothesis, we applied a broken PL 
model to the PSPC and ACIS spectra and obtained an accept¬ 
able fit (x^ = 1.31) with nH ,20 = 0 - 3 !! 3 j 2 and photon indices 
7 i = 2.0±0.2 and 72 = 3.6 ±0.1, below and above the break 
energy Eb= 1.1 ±0.1 keV, respectively. However, the EUVE 
DSI flux calculated with the best-fit parameters exceeds the 
observed DSI flux by about 3cr (if one adopts a 15% system¬ 
atic uncertainty in the DSI flux). Including the DSI data in the 
combined analysis shows that the broken PL model can be ac¬ 
cepted at the same low-confidence level as in the case of the 
PL plus single-temperature polar cap model fit, and it requires 
nH ,20 > 0.5. 

®See: 


Broken PL spectra have been observed from young ordinary 
pulsars, but with a much higher break energy, in the GeV range. 
However, no breaks have been observed in the spectra of other 
millisecond pulsars, and the (single) PL slopes measured in the 
three millisecond pulsars with certainly nonthermal radiation 
are more gradual than the slopes of the broken PL components 
in PSR 10437^715. Therefore, although we cannot completely 
reject the broken PL interpretation, we consider it less plausible 
than the power-law plus thermal polar-cap model. 


4. HRC OBSERVATIONS AND TIMING ANALYSIS 

PSR 10437^715 was observed with the imaging array of the 
High Resolution Camera (HRC-I; Murray et al. 1997) on 2000 
February 13, with a 18.9 ks effective exposure. The pulsar was 
detected with a source count rate of 174 ±4 counts ks“' (es¬ 
timated in a 3"-radius circle centered at the pulsar’s position). 
The image of the pulsar and its surroundings is shown in Figure 
3. The surface brightness of the background in the pulsar vicin¬ 
ity, 0.015 counts ks“' arcsec“^ is so low that its contribution is 
negligible in the source aperture, and the count rate error is de¬ 
termined solely by the source photon statistics. Figure 4 shows 
the extracted radial distribution of the (energy-integrated) encir¬ 
cled source count fraction. The distribution is well consistent 
with that expected for a point source. 

Our attempts to detect pulsations in these data failed. The 
reason of this failure was found later — a wiring error in the 
HRC causes the time of an event to be associated not with that 
event, but with the following event at the HRC front end, which 
most likely is a background particle event rejected by the anti- 
coincidence shield and not telemetered to the ground (see POG, 
§7.10; Tennant et al. 2001). As a result, most of the true event 
times are lost, and the HRC timing accuracy degrades from the 
planned 16 ps to the mean time between events at the HRC 
front end, typically about 4 ms. 

To circumvent this problem, the HRC team worked out a 
special operating mode which employs the HRC-S instrument. 
This mode uses only the central segment of the HRC-S, with 
two outer segments disabled, which lowers the overall counting 
rate below the telemetry saturation limit, even with all on-board 
vetoes disabled. As a result, all events are telemetered and can 
be assigned correct times by shifting the time-tags in ground 
processing. 

The pulsar was re-observed in this new mode on 2000 Oc¬ 
tober 6 (19.6 ks effective exposure). We measured the source 
count rate of 152 ±6 counts ks“' using the same aperture as 
for the HRC-I data. To perform the timing analysis, we first 
extracted the “arrival times” (as given in Level 2 data file) of 
2,795 events from a l"-radius circle centered at the pulsar (93% 
of the total number of source events). Then, for each of the 
selected events, we found the time of the following event in 
the Level 1 data file (which, in addition to photon events, con¬ 
tains time-tagged particle events) and assigned this later time 
to the selected source event. These corrected times were con¬ 
verted to equivalent arrival times at the solar system barycenter 
(SSB) in the barycentric dynamical time (TDB) system with the 
axBary tool® using the solar system ephemeris DE200 (Stan- 
dish 1982). 

To correct the TDB arrival times for the effect of the orbital 
Doppler shift (Taylor & Weisberg 1989) and extract the light 


http://asc.harvard.edu/ciao2,2/test/ahelp/axbary.html 
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curve, we used the pulsar and binary orbit ephemerides given in 
Table 1. The values of statistics (Buccheri et al. 1983) calcu¬ 
lated with the X-ray arrival times at the radio ephemeris param¬ 
eters are Z\ = 217, Z| = 258, Z| = 259. The difference between 
the radio period and that corresponding to the Z| maximum, 
6P = 3x 10“*^ s (or SF = dP/P^ = 1 x 10“^ Hz), is negligible for 
the timing analysis of the X-ray data collected in a 20 ks time 
span. Note that SP is much smaller than the 1 a uncertainty 
(given by the Bayesian method — see Zavlin et al. [2000] for 
details), 1 x 10“^“ s, of the period found in the HRC-S data. The 
small value of SP and the large Z^ values prove unambiguously 
that the new HRC-S timing mode is fully capable of detecting 
millisecond pulsations. (The largest Z^ values achieved in the 
ROSAT observations of this pulsar are 117 and 125 for n = I 
and n = 2, respectively). 

The upper panel of Figure 5 presents the extracted HRC- 
S (energy-integrated) pulse profile of PSR J0437^715. Us¬ 
ing the radio ephemeris given in Table 1 and radio observa¬ 
tions at frequencies close to 1400 MHz, made using the Parkes 
Telescope and the Caltech correlator (cf. Sandhu et al. 1997), 
which were contemporaneous with the X-ray observations, we 
were able to find the relation between the radio and X-ray 
phases — (/) = 0 in this figure corresponds to the maximum 
of the radio peak observed at 1420 MHz, dedispersed to infi¬ 
nite frequency using the DM in Table 1. The epoch of (/> = 0 
is 51823.7311148130996 MJD [UTC at geocenter]. Using 
tempo'®, we determined the phase offsets of the X-ray peaks 
for five sub-intervals within the 20-ks exposure. The corre¬ 
sponding time offsets vary from -6 /rs to h- 37 /rs, with a mean 
value of h-16 ^s (0.3% of the pulsar period), equal to the in¬ 
trinsic HRC time resolution. Therefore, we conclude that the 
phases of the X-ray and radio pulses are virtually the same — 
their difference does not exceed the errors of the Chandra abso¬ 
lute time". It should be noted that this is the first determination 
of the absolute phase of X-ray pulsations of PSR J0437-4715 
— errors in absolute timing of the ROSAT data were at a 1-2 
ms level (see HMM96 and BT99 for discussion), preventing the 
absolute phase determination. 

For comparison, we also show in Figure 5 the pulse profiles 
obtained from the ROSAT PSPC and HRI observations of 1994 
(with phases chosen to approximately co-align the ROSAT 
peaks with Chandra peak). The pulsed fraction of the HRC-S 
light curve is somewhat larger than those in the ROSAT data. 
Since the HRC-S is more sensitive at higher energies than the 
ROSAT detectors, this may indicate that /p increases with E, 
as predicted by the hydrogen polar cap model. In addition, the 
pulse in the HRC-S light curve looks slightly narrower than the 
PSPC and HRI pulses. This also could be explained by the 
properties of the thermal radiation emitted from a hydrogen at¬ 
mosphere — the radiation intensities are more beamed at higher 
E (Zavlin et al. 1996). 


5. DISCUSSION 


The Chandra observations of PSR J0437^715 have allowed 
us to perform the spatial, timing and spectral analyses of the 
new data collected with high angular and spectral resolution in 
an extended energy range. 

The principal new result of the ACIS observation is the mea¬ 
surement of the pulsar’s X-ray spectrum at higher energies, up 
to 7 keV. Among the formally acceptable fits of the combined 
ROSAT PSPC and Chandra ACIS spectra, only the broken PL 
model corresponds to a purely nonthermal pulsar radiation in 
the X-ray band. However, the break energy of about 1.1 keV 
is 5-6 orders of magnitude lower than those observed in other 
pulsars. Such a difference looks too large to be explained by 
a lower magnetic field in the radiating region — it would also 
require a lower break energy in the spectrum of radiating elec¬ 
trons/positrons, compared to ordinary pulsars. Moreover, ex¬ 
trapolation of the broken PL model to the optical B and V 
bands (assuming the same slope as in the X-ray range below 

1.1 keV) predicts optical magnitudes, = 19.6 ± 1.4 and my = 

19.1 ±1.4 (for extinction coefficients Ab = 0.3, Ay = 0.2 — 
Danziger, Baade, & della Valle 1993), much brighter than those 
detected from the white dwarf companion, = 22.1 ± 0.1 and 
my = 20.9 ±0.1 (Bailyn 1993; Danziger et al. 1993; Bell et 
al. 1993). Thus, we do not consider the broken PL model to be 
a plausible interpretation. 

The alternative model involves two components of different 
origin — thermal and nonthermal. The nonthermal component 
originates in the pulsar magnetosphere'^, whereas the thermal 
component is emitted from hot polar caps on the NS surface. 
Depending on assumption about the polar cap temperature dis¬ 
tribution, one gets different relative contributions from these 
two components. In the model with uniformly heated polar 
caps, the nonthermal component described as a PL of photon 
index 7 = 2.7-2.9 provides about 80% in the X-ray flux and 
dominates at energies below 0.6 keV and above 2.7 keV. How¬ 
ever, this model encounters the same problems in the EUV and 
optical bands as the broken PL (the steeper slope of this PL 
component predicts even higher fluxes in the B and V bands). 

If we assume a more plausible polar cap model with temper¬ 
ature decreasing outwards from the cap center, then the ther¬ 
mal component becomes dominant between 0.06 keV and 2.5 
keV, providing some 75% of the X-ray flux, while the PL com¬ 
ponent of 7 = 1.6-2.5 dominates outside this band. In ad¬ 
dition to a more realistic temperature distribution, the latter 
model is well consistent with the EUVE data and yields esti¬ 
mates on the hydrogen column density in agreement with the 
indirect measurements. As this PL component is fainter than 
in the two other models, its extension to the optical falls be¬ 
low the observed radiation of the white dwarf companion for 
photon indices 7 < 1.9. This allows us to predict that the PL 
component should be observable in the UV (particularly, far- 
UV) range where it is brighter than the Wien tail of the white 
dwarf spectrum, assuming there is no turnover of the nonther¬ 
mal spectrum between the UV and soft-X-ray energies. On the 
other side of the X-ray band, extrapolation of the PL compo¬ 
nent to the gamma-ray energies above 100 MeV predicts a pho- 


‘OTEMPO k a program fnr the analysis nf pnkar timing Hafa maintainpH anH HktrihnfpH hy Prinrpton University and the Australia Telescope National Facility. It 
is available at http : / /www. atnf. csiro . au/re search/pulsar/timing/ tempo /[ 

^^The accuracy of the Chandra absolute time was estimated by Tennant et al. (2001) from the comparison of the Chandra HRC-I and RXTE timing of the Crab 
pulsar: fchandra-?RXTE = “200 ib 100 /TS. 

^^Conventionally defined as the ratio of the number of counts above the minimum of the light curve to the total number of counts. 

^^Another source of unresolvable nonthermal radiation could be the shocked pulsar wind near the white dwarf companion (e.g., Arons & Tavani 1993). However, 
at the distance of ap = 1.5 x 10^^ cm from the pulsar (van Straten et al. 2001), the companion intercepts only a fraction 6 x 10“^ of the wind (assuming the wind is 
approximately isotropic), too small to explain the observed nonthermal luminosity. 
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ton flux / < 2 X 10“* s“^ cm“^ (7 > 1 -6), below the upper limit, 
/ < 1.5 X 10“^ s“^ cm“^, obtained from the CGRO EGRET ob¬ 
servations (Eierro et al. 1995). 

We emphasize that these models require the thermal radia¬ 
tion to be emitted from hydrogen (or helium) NS atmosphere. 
The high spectral resolution of the ACIS data rules out an at¬ 
mosphere comprised of heavier chemical elements. 

The HRC-S observation of PSR J0437^715 has demon¬ 
strated the Chandra timing capability at a millisecond level. 
The HRC-S pulse profile looks narrower, and the pulsed frac¬ 
tion is somewhat higher, than those obtained in the earlier 
ROSAT and EUVE observations at lower energies, which 
could be explained by the properties of the thermal radiation 
from polar caps covered with a hydrogen or helium atmosphere. 
On the other hand, the shape of the profile is clearly asym¬ 
metric, with a longer rise and faster decay, which cannot be 
explained by a simple axisymmetric temperature distribution. 
Relativistic effects (particularly, the Doppler boost) should lead 
to a different asymmetry — a faster rise and longer trail (Braje 
& Romani 2000; Eord 2000). The analysis of HRC-S data 
demonstrates, for the first time, that the phase of the X-ray pulse 
virtually coincides with that of the radio pulse. If, as we sug¬ 
gest, the main contribution to the HRC-S band is due to the 
thermal polar cap radiation, and if the pulsar radio beam is di¬ 
rected along the magnetic axis, then the radio emission must be 
generated close to the NS surface — e.g., the time difference 
of <0.1 ms between the X-ray and radio phases corresponds 
to a distance of < 30 km, much smaller than the light cylin¬ 
der radius, R\c = 275 km. Alternatively, if the radio emission 
is generated at a higher altitude, the combination of field-line 
sweepback and aberration must contrive to cancel the radial 
travel-time difference. 

The Chandra observations show no sign of an X-ray PWN 
that could accompany the bow-shock revealed by the H^ ob¬ 
servations. Three-sigma upper limits on the PWN bright¬ 
ness (in counts arcsec”^) can be estimated as 3(h/A)^/^, where 
b is the background surface brightness (b = 0.51 and 0.28 
counts arcsec“^ for the ACIS-S and HRC-I images, respec¬ 
tively), and A is the PWN area (we will scale it as A = 1000/a 
arcsec^, assuming that a typical transverse size of the PWN 
somewhat exceeds the stand-off distance, 10", of the bow 
shock). Eor a power-law PWN spectrum with a photon index 
7=1.5-2 (similar to those observed from other PWNe), these 
upper limits correspond to the PWN intensities / ^wn < (13- 

1.8) X 10-'V;‘/" and /.p^n < (3.6-5.7) x lO-'V/'^" erg cm-^ 
s“' arcsec”^, for the ACIS-S and HRC-I, respectively, in the 
0.1-10 keV range. The corresponding upper limits on the PWN 
X-ray luminosity, L^^pwn ~ 47rc/^A4^pwn are much smaller than 
the rotational energy loss rate, E = 3.8 x 10^^ erg s“' — e.g., 
Ex,pwn < (3.0^.2) X 10^*^/J^^ erg s“' for the more sensitive 
ACIS-S limit. 

The low upper limits on the PWN luminosity in X-rays 
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can be simply explained by a low magnetic field in the PWN 
region, expected for a particle-dominated pulsar wind. The 
shock in the relativistic pulsar wind should be located just in¬ 
terior to the observed Hq, bow shock (Arons & Tavani 1993). 
When the wind electrons pass through the shock, their direc¬ 
tions of motion become “randomized”, and their synchrotron 
radiation may result in an X-ray nebula, provided the electron 
energies and the magnetic field are high enough in the post¬ 
shock region. The pre-shock magnetic field can be estimated as 
Bi = [£/(/or2c)]'/2[a/(l + a)]'/2= 18/-'/'[a/(l + a)]'/2 ^iG, 
where « 2 x lO'® cm is the stand-off distance corresponding 
to 10" atd= 140 pc, /o = AH/(47r) < 1 is the collimation fac¬ 
tor of the wind, and the “magnetization parameter” a is the ratio 
of the Poynting flux to the kinetic energy flux. The maximum 

- 1/2 

value of the post-shock magnetic field, B ~ B 1 ~ 18 pG, is 
obtained for cr ^ 1. However, according to Kennel & Coronity 
(1984; KC84 hereafter), a significant fraction of the total energy 
flux upstream can be converted into (observable) synchrotron 
luminosity downstream only if a <0.1 (e.g., these authors esti¬ 
mate a w 0.003 for the Crab pulsar). Eor cr <C 1, the post-shock 
magnetic field is B ^ 3(l-4cr)Bi ~ 53/f2*^^cr'/^(l-4.5 ct) /xG 
(e.g., from eqs. [4.8] and [4.11] of KC84, B increases from 
3 /iG to 12 /iG when a increases from 0.003 to 0.1, at /o = 1). 
Such low magnetic fields in the shocked wind strongly limit 
the maximum energy, mgC^Emax, of radiating electrons and, 
consequently, the maximum frequency t'max of the synchrotron 
radiation. Since the Larmor radius of most energetic elec¬ 
trons, rL= 1.7 X cm, cannot exceed the shock 

radius substantially, we obtain Emax < lO^/^B-s, ht'max ~ 
{heB/4TTmec)Tl^^^ < 0.6//B_g keV, where B-s = B/(10/xG), 

fs = ri/rs ^ If, for instance, B < fiG (i.e., a < 

-4/3 

O.Ol fs in the KC84 model), the synchrotron emission at the 
bow shock is not expected to be seen in X-rays. 

Despite the superior quality of the Chandra data, which have 
allowed us to detect the hard tail of the pulsar’s spectrum and 
pinpoint the absolute phase of the X-ray pulse, there still re¬ 
main some open problems. Although our analysis strongly fa¬ 
vors the thermalH-nonthermal interpretation, it is still unclear 
which of the two components dominates in the X-ray radia¬ 
tion of PSR 10437^715. To establish the relative contributions 
of these components, energy-resolved timing and time-resolved 
spectral analysis are needed, which, hopefully, will be possible 
with the forthcoming XMM-Newton data. 

We thank Leisa Townsley for the advise on the ACIS data 
reduction, Allyn Tennant for the discussion of the HRC tim¬ 
ing issues, and Willem van Straten for providing the timing 
ephemeris given in Table 1. We also thank Vadim Burwitz for 
helping in preparation of the color image. This work was par¬ 
tially supported by NASA through grants NAG5-7017, NAG5- 
10865, and SAG GOO-1126X. 
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Table 1 


Ephemeris parameters for the binary pulsar 30437^715“ 


Parameter 

Value 

R.A. (J2000). 

04''37”15?78651377 

Dec. (J2000). 

-47°15' 8"46158135 

(mas/yr).. 

. 121.4565 

fis (mas/yr). 

-71.4494 

Parallax (mas). 

. 6.9830 

F (Hz). 

173.68794899909872242 

F (10“'^ s-2). 

-1.728325769909 

F (10-^° s-3). 

132.1214177 

Epoch of period and position (MJD). 

51,194.0 

DM (cm“^ pc). 

2.64690 

Binary period, P\, (days). 

5.741044490290 

x = flp sin / (It-s)?.. 

3.366693427 

Eccentricity, e . 

0.000019193 

Longitude of periastron, oj (deg). 

1.848063 

Epoch of periastron (MJD). 

. 51,194.634254236 

i(10“i2 It-s s-'). 

0.078270 

Pb (10“'2 s s-'). 

, 3.397 


Based on radio observations made using the Parkes 64-m telescope of the Australia Telescope National Facility (van Straten et al. 2001) 
Projected semi-major axis of the pulsar’s orbit 
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Fig. 3.— 5'.2 x 6.'3 region of the HRC-I image. The bright source in the upper-left corner is an AGN (Halpern & Marshall 1996). 
The inset in the upper-right corner shows a zoomed 63" x 63" field around the pulsar. The dashed curve marks the leading edge of 
the Hq, bow-shock. The arrow shows the direction of the pulsar’s proper motion, with the transverse velocity v±. 
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Rotational phase 0 

Fig. 5.— Pulse profiles of PSR J0437^715 from the Chandra, ROSAT, and radio observations. The smooth light curve in the upper 
panel is obtained by averaging of the 20-bin profile over the starting phase. The zero phase corresponds to the peak of the radio pulse 
shown in the lower panel which is the summed profile from the 1400 MHz radio observations made with the Caltech correlator at 
Parkes. The absolute phases of the ROSAT pulses were not measured; we chose the phases to approximately co-align these pulses 
with the Chandra pulse. 




